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Certain real parameters of a Hamiltonian, when continued to complex values, can give rise to 
singular points called exceptional points (JfP’s), where two or more eigenvalues coincide and the 
complexified Hamiltonian becomes non-diagonalizable. We show that for a generic d-dimensional 
topological superconductor/superfluid with a chiral symmetry, one can find EP’s associated with 
the chiral zero energy Majorana fermions bound to a topological defect/edge. Exploiting the chiral 
symmetry, we propose a formula for counting the number (n) of such chiral zero modes. We also 
establish the connection of these solutions to the Majorana fermion wavefunctions in the position 
space. The imaginary parts of these momenta are related to the exponential decay of the wavefunc¬ 
tions localized at the defect/edge, and hence their change of sign at a topological phase transition 
point signals the appearance or disappearance of a chiral Majorana zero mode. Our analysis thus 
explains why topological invariants like the winding number, defined for the corresponding Hamilto¬ 
nian in the momentum space for a defectless system with periodic boundary conditions, captures the 
number of admissible Majorana fermion solutions for the position space Hamiltonian with defect(s). 
Finally, we conclude that EP’s cannot be associated with the Majorana fermion wavefunctions for 
systems with no chiral symmetry, though one can use our formula for counting n, using complex k 
solutions where the determinant of the corresponding BdG Hamiltonian vanishes. 

PACS numbers: 73.20.-r, 74.78.Na, 03.65.Vf 


Introduction - The Hamiltonian operator can con¬ 
tain certain real parameters, which on being continued 
to complex values, give rise to singular points where the 
operator becomes non-diagonalizable. These are called 
exceptional points {EP’s), at which two or more of the 
eigenvalues coalesce and the norm of at least one eigen¬ 
vector of the complexified Hamiltonian vanishes 
The concept of EP’s is similar to that of a degeneracy 
point, but with the important difference that all the en¬ 
ergy eigenvectors cannot be made mutually orthogonal. 
In previous works, EP’s have been used [EHS] to describe 
topological phases of matter associated with zero energy 
Majorana bound states (MBSs) in one-dimensional (Id) 
topological superconductors/superfluids. 

The first-quantized Hamitonians describing fully 
gapped noninteracting topological insulators and super¬ 
conductors in d-dimensions can be classified into ten 
symmetry classes m I14j in terms of nonspatial sym¬ 
metries, i.e., symmetries that act locally in the position 
space, namely time-reversal symmetry (TRS), particle- 
hole symmetry (PHS), and chiral symmetry. Recently, 
it has been realized [THtiST] that the complete classifica¬ 
tion should include topological states protected by crys¬ 
talline symmetries (such as mirror reflections and rota¬ 
tions), which are spatial symmetries acting nonlocally in 
the position space. 

A superconductor is described by a Bogoliubov de 
Gennes (BdG) Hamiltonian (idsdc)) which has an exact 
PHS, in addition to having the structure of an ordinary 
Bloch Hamiltonian. The PHS of H^dc implies 

CH^do{k)C-^ = -H*^^{-k), C = (1) 


where C is an antiunitary operator, so that the energy 
eigenvalues ±E{k) always appear in pairs for each k- 
value. In this paper, we will consider the case = 1 
for which idsdc can be categorized into three non-trivial 
classes: BDI, Dili and D. 

Topologically protected gapless modes can occur at a 
topological defect or surface of a topological supercon¬ 
ductor/insulator. Due to the PHS, the boundary ex¬ 
citations of a topological superconductor are Majorana 
fermions, such that the creation operator (yj^^p) is equiv¬ 
alent to the annihilation operator {je-o)- Since they 
obey non-Abelian statistics, they might have potential 
applications in designing fault-tolerant quantum comput¬ 
ers 

In an earlier work [T^, a formula to count the Ma¬ 
jorana zero modes in Id chiral topological superconduc¬ 
tors/superfluids based on the notion of EP’s has been 
proposed. This formula serves as an alternative to the 
familiar winding number (W) [T31II111M1- In this paper, 
we will show that we can apply the EP formalism to the 
zero energy Majorana excitations bound to defects for 
generic d-dimensional Hamiltonians anticommuting with 
a chiral symmetry operator. We will also generalize the 
formula to count the number (n) of Majorana fermions 
in the presence or absence of the chiral symmetry. 

EP solutions for systems with a chiral sym¬ 
metry - We consider a topological defect embedded 
in (or at the boundary of) a d-dimensional topological 
superconductor. Let m be the dimensions of the de¬ 
fect, parametrized by the Cartesian coordinates r± = 
{ri,...,rd-m) and ry = {rd-m+i, ■ ■ ■ ,rd), located at 
r± — 0. Let k± = k±fl = {ki,... ,kd-m) and fey = 
{kd-m+i, ■ ■ ■ ,kd) be the corresponding conjugate mo- 
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menta, where fc_L = |fcx| and is the unit vector when 
written in spherical coordinates. If there is a chiral sym¬ 
metry operator O which anticommutes with the Hamil¬ 
tonian, the latter takes the form 


{ 0 A{-^V) \ 


( 2 ) 


transition point, the position space wavefunction of the 
corresponding Majorana fermion changes from exponen¬ 
tially decaying to exponentially diverging or vice versa. 
If the former happens, the Majorana fermion ceases to 
exist. A new Majorana zero mode appears in the lat¬ 
ter case. The count (n) for the Majorana fermions for a 
defect is thus captured by the function 


in the eigenbasis of O (where it is diagonal). Here A 
is an N X N matrix operator in the position space. The 
solutions for the Majorana zero modes are given by "0+ = 
(m+, 0)^ and V'- = (0,u_)^ (with chirality -|-1 and —1 
respectively, since 0'ip± = ±ip±, satisfying: 

A{-iV)u-=0, M^(-iV)u+ = 0. (3) 

Assuming the dependence U- ^ exp (—z°‘r°‘) in the bulk, 
we find that z“’s should satisfy 

det {A{iz)) = 0 . (4) 

Now let us find out the EP solutions for the Hamilto¬ 
nian 

T)- <"> 

in the momentum space, for the corresponding system 
with no defect. On analytically continuing the magnitude 
= k = \k\to the complex fcx-plane, at least one of the 
eigenvectors of iJchirai(fc) collapses to zero norm where 

det {A{k)) =0, or det (M^(fc)) =0. (6) 

These points are associated with the solutions of EP’s 
for complex fcx-values where two or more energy levels 
coalesce. Furthermore, these coalescing eigenvalues have 
zero magnitude as det(M(fc)) = 0 (or det (M^(fc)) = 0) 
also implies det (i4„hirai(fc)) = 0. Hehirai(fc) becomes non- 
diagonalizable, as in the complex A:x-plane, det (M(fe)) = 
0 ^ det (M^(A;)) = 0 (or vice versa). However, at the 
physical phase transition points, the imaginary parts 
of one or more solutions vanish, and det(M(A;)) = 
det(A^(fe)) = 0 for those solutions, making i4„hirai(fe) 
once again diagonalizable and marking the disappearance 
of the corresponding AP’s. We find that both k± and 
iz (where z = \z\) are obtained as the solutions of the 
same equation (Eq. immediately indicating the cor¬ 
respondence ik± —z. Each EP solution corresponds 
to a Majorana fermion of a definite chirality. 

For a generic Pchirai, let k\ and (j = 1,..., Q) 
be the two sets of EP solutions for det {A{k)) = 0 and 
det(M^(fc)) = 0 respectively, related by {Im(fc')^)} = 
— {Im(fc;^)}, after k± has been analytically continued to 
the complex plane. At a topological phase transition 
point, one or more of the Im go through zero. 

When Im {k\^g^ changes sign at a topological phase 


1 ^ 
i=i 

-S5n{ Im ({A°},fc[|,n°)) } , 

(7) 

where ({Ai}, fey, FI) are the parameters appearing in the 
expressions for and ^{A°}, fc{j, 17°^ are their values 

at any point in the non-topological phase. If Al^(fc) = 

A"^{—k) holds, then the two sets of EP's are related by 
{k\} = —{kPg}, one set corresponding to the the solu¬ 
tions obtained from one of the two off-diagonal blocks. 

In such cases, the pairs of the Majorana fermion wave- 
functions are of opposite chiralities. 

This connection between the complex /c^-space solu¬ 
tions where two or more zero energy eigenvalues coa¬ 
lesce, with the Majorana fermion solutions of the po¬ 
sition space Hamiltonian with a topological defect, helps 
us understand why the topological invariants like W 
(winding number [241) /({^i}) ^|| > ^)) defined for 

systems with a chiral symmetry, are related to n. W 
and /({Ai}, fc||, J7) are defined in terms of det {A{k)) for 
TIchirai written in the momentum space, in a given topo¬ 
logical phase. On the other hand, Majorana fermion 
solutions are calculated for the corresponding Hamilto¬ 
nian at the location of the defect. Along an EP so¬ 
lution, the complexified Hamiltonian has vanishing de¬ 
terminant, indicating the presence of two or more coa¬ 
lescing zero eigenvalues. The actual physical Hamilto¬ 
nian, however, has zero determinant only for a real k±- 
value, at which one or more AP’s collapse, allowing the 
Hamiltonian to be diagonalizable at that point. Observ¬ 
ing the correspondence exp^ifc^y^^ exp(—z^-rx), 

such that ’’-l) is a Majorana fermion so¬ 

lution, explains why a Majorana zero energy state exists 
throughout a topological phase, with the topological in¬ 
variants /({Ai}, fc||, f7) and W capturing the number of 
such modes, and changing their values only on crossing 
a phase transition point. This also provides an intuitive 
understanding of the bulk-edge correspondence. 

EP solutions for systems without chiral sym¬ 
metry - A Hamiltonian Hg for a system without any 
chiral symmetry can be constructed from Hchirai as [25j 

Hg = cos 9 P.hirai + sin9 0 , 
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This clearly breaks the chiral symmetry unless 9 = 0. 
The EP’s are given by det {A{k)) = 0 or det (^'^(fe)) = 
0, where two levels coalesce for a complex value of k±. 
However, we immediately observe that these EP^s do 
not correspond to zero energy modes for the complex¬ 
ified In the momentum space, at the points 

of vanishing energy eigenvalues (and hence vanishing 
det (i?D(A;))), Hoik) still remains perfectly diagonaliz- 
able even for complex /c_l- values. 

The position space Majorana fermion solutions thus 
cannot be captured by the EP solutions in the complex 
A:_l- space, as the latter are related to vanishing eigenvec¬ 
tors for the complexified Hamiltonian where it becomes 
non-diagonalizable. However, we can still encode the 
count of the non-chiral Majorana zero modes by the func¬ 
tion defined in Eq. but with being the two sets 
of complex-valued solutions for det (Hoik)) = 0, obeying 
{lm(k\)} = -{Im(A:^)}. 

Point Defects - In this section, we consider point 
defects for a d-dimensional BdG Hamiltonian with the 
PHS operator C squaring to -I-1. Such a defect can ap¬ 
pear at the edges of a Id chain, as vortex solutions in a 
2d system, or as hedgehog configurations in a 3d bulk. 
There can be a Majorana zero energy state bound to 
such a topological defect, depending on the values of the 
parameters. 

Point Defects in class BDI - The BDI symmetry 
class is associated with the existence of topologically pro¬ 
tected chiral Majorana fermions and characterized by an 
integer (Z) topological invariant. 

We consider a specific model [55] for such point defects 
described by a BdG-Dirac type Hamiltonian: 

= -ir“ , (9) 

where a, (3 = I,...,d, (jip is a component of the gen¬ 
eralized d-component order parameter (p, and the 2d- 
dimensional gamma matrices obey the anticommutation 
relations {r^,r‘^} = 2 5^'^. Clearly, r“ and T^ are the 
sets of gamma matrices associated with the kinetic term 
and the order parameter part of idchira,i- 

The point defect can be modelled by choosing cp of the 
form: 


Pair) = A{r)r°‘, A(r) 


0 for r < i? , 
Ao for r > i?, 


( 10 ) 


where Aq is a constant. The chirality operator O can be 
identified with 


p2d-Hl ^ (11) 

which anticommutes with all the T'^’s. 

One can choose the gamma matrices as follows |40| : 

r“ = 7“ 0 1 (g) (Ti, r‘^+“ = I (g 7“ 0 (72, 

= 1010 (73. (12) 


The normalizable solutions |26j are given by 




= (7+ ( Q j exp [(- 1 ) y dr'A(r')J for (-l)^Ao < 0, 




= C_ ( j, ) exp 


(-1)^+1/ dr'A(r') for (-l)^Ao >0, 

Jo J 

(13) 


where C± are constants, and £ and T have been defined 
in Eq. (16). Hence the admissible solutions decay as 


exp(—I AqI rj_) in the bulk, where i z = i Aq can be iden¬ 
tified with the EP solution with complex k±. 

In Id, numerous lattice models have been studied 
which can support one [55] or multiple MBSs |57H55] at 
one edge of an open chain. Eor such models, one can 
show that the chiral MBS solutions at an edge, obtained 
by the transfer matrix approach formalism [5511571 - 1551155] 
in the lattice space, have the same relation with the EP 
solutions obtained in the complex fc-space. 

Point Defects in class D - In the symmetry class D, 
the TRS is broken, thus breaking the chiral symmetry, 
and characterized by a Z 2 topological invariant |25j , with 
n allowed to be 0 or I. A topologically protected non- 
chiral MBS bound to a point defect can exist. In the 
absence of O, the conclusions iov Ho in Eq. ([^ will apply. 

Point defects in class Dili - A point defect in class 
Dili can support a Majorana Kramers pair (MKP) cor¬ 
responding to doubly degenerate Majorana zero modes, 
and is characterized by a Z 2 topological invariant. A chi¬ 
ral symmetry operator O can be defined such that the 
Hamiltonian in class Dili can be brought to the block 
off-diagonal form (Eq. just like for the class BDI. 
Eurthermore, the Dili class with a mirror symmetry is 
equivalent to the BDI class with an additional (pseudo) 
TRS, and hence the edge of a Id Hamiltonian with a 
mirror line can be characterized by an integer (Z) mir¬ 
ror winding number [IMIj related to the existence of 
multiple MKPs. Hence the analysis and conclusions for 
ddchirai in Eq. ([^ will be applicable for the Dili class. Due 
the connection between EP's and the Majorana fermion 
wavefunctions, we can again use Eq. 0 to count the 
number of MKPs in a given topological phase in such 
a Id system. Each Majorana zero mode is of a definite 
chirality with respect to O. 

Line Defects - In this section, we consider line de¬ 
fects for a d-dimensional BdG Hamiltonians with = 1. 
Majorana zero modes bound to such a defect can appear 
at the edges of a 2d system, or along a vortex line in a 
3d system [51] . 

Whenever H^dc cannot be written in a block-off diago¬ 
nal form in the absence of a chiral symmetry, the system 
will conform with the discussion for Ho vn Eq. ([^. 

For a system in class BDI or class Dili, the chi¬ 
ral symmetry operator O exists, so that the Hamil- 
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tonian -ffchirai takes the block off-diagonal form shown 
in Eq. ([^. The allowed solutions will be of the form 
^ exp (i fc|| ry) exp (—z |r_L|). Clearly, i z will correspond 
to the EP solutions in the complex fcj_-plane satisfying 
Eq.(|^. Using these solutions, /({Ai}, /cy, fi) in Eq. Q 
will give the count of the chiral Majorana fermions. 

Conclusion - We have established the relation of 
the EP solutions for complexified momenta to the Ma¬ 
jorana fermion wavefunctions bound to a topological de¬ 
fect, for a generic d-dimensional topological superconduc¬ 
tor having a chiral symmetry operator O. This connec¬ 
tion explains why topological invariants such as the wind¬ 
ing number, defined for the corresponding BdG Hamilto¬ 
nian without the defect (or edge) written in the momen¬ 
tum space, capture the number of admissible Majorana 
zero mode solutions for the position space Hamiltonian 
describing the defect. Each of these solutions are of a 
definite chirality (with respect to O). We have also pro¬ 
posed a formula for counting the number of Majorana 
zero modes at each defect/edge, based on these EP so¬ 
lutions. We have shown that such EP solutions cannot 
exist for systems without any chiral symmetry (e.g. class 
D). Despite such solutions not existing in the absence 
of O, we have, nevertheless, established an expression 
for calculating the non-chiral Majorana fermions at a de¬ 
fect from the solutions of det {El^^ci^)) = 0 in the com¬ 
plexified momentum space. Though we have explicitly 
discussed only point and line defects, the arguments pre¬ 
sented will hold for any m-dimensional defect embedded 
in a d-dimensional superconductor (with d > m), or for 
the (d— l)-dimensional boundary of the bulk system. A 
variety of Id and 2d models has been studied in |35j . ex¬ 
plicitly demonstrating how our counting formula works. 
A proof of the counting formula in Eq. Q has also been 
been provided in this follow-up work. Einally, it will be 
interesting to explore the connection between our count¬ 
ing formula with those proposed in earlier works [55H55] . 
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